Diffusion and dispersion of passive tracers: Navier-Stokes versus MHD turbulence 
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A comparison of turbulent diffusion and pair-dispersion in homogeneous, macroscopically isotropic 
Navier-Stokes (NS) and nonhelical magnetohydrodynamic (MHD) turbulence based on high- 
resolution direct numerical simulations is presented. Significant differences between MHD and NS 
systems are observed in the pair-dispersion properties, in particular a strong reduction of the sepa- 
ration velocity in MHD turbulence as compared to the NS case. It is shown that in MHD turbulence 
the average pair-dispersion is slowed down for Td < t < lOYd, Td being the Kolmogorov time, due 
to the alignment of the relative Lagrangian tracer velocity with the local magnetic field. Signifi- 
cant differences in turbulent single-particle diffusion in NS and MHD turbulence are not detected. 
The fluid particle trajectories in the vicinity of the smallest dissipative structures are found to be 
characterisically different although these comparably rare events have a negligible influence on the 
statistics investigated in this work. 

PACS numbers: 47.10-g;47.27.-i;52.30.Cv 

The diffusive effect of turbulence on contaminants passively advected by the flow is of great practical and fundamen- 
tal interest. While the study of passive scaiars [l| usually reverts to the Eulerian description of the flow, the Lagrangian 
point of view has proven to be very fruitful regarding investigations of turbulent diffusion and pair-dispersion 0, 0] 
as well as for the fundamental understanding of turbulence [4|. The three-dimensional dynamics of passive tracers 
in neutral fluids has been subject of various experimental (for recent works see e.g. 0, H, 0]) an d numerical, e.g. 



lfj,lll|, investigations. Related problems regarding the turbulent diffusion of magnetic fields and the influence 
of turbulent magnetic fields on particle diffusion have been studied extensively in space and astrophysics, see e.g. 
12 , 15 . Im 17|. as well as in the context of magnetically confined nuclear-fusion plasmas, see for example 

This Letter reports a first effort to identify differences in the diffusion and dispersion properties of turbulent flows 
in electrically conducting and in neutral media. To this end the dynamics of fluid particles is studied via high- 
resolution direct numerical simulation of passive tracers immersed in fluids that are described by the incompressible 
magnetohydrodynamic (MHD) and the Navier-Stokes (NS) approximation. 

Using the vorticity, oj = V x v, and a uniform mass density, po = 1, the non-dimensional incompressible MHD 
equations are given by 

d t us = Vx[t)xu-h(Vxfi)] + fiAu) (1) 
d t b = V x (v x b) + -qAb (2) 
V-v = V-6 = 0. (3) 

The dimensionless molecular diffusivities of momentum and magnetic field are denoted by /i and 77, respectively. The 
magnetic field b is given in Alfven-speed units. The Navier-Stokes equations which govern the motion of an electrically 
neutral fluid are obtained by setting b = in Eqs. (HJ)-©. 

The MHD/Navier-Stokes equations are solved by a standard pseudospectral method in a triply periodic cube of 
linear extent 2tt. The velocity field at the position X of a tracer particle is computed via tricubic polynomial 
interpolation and then used to move the particle according to 

X(t) = V(t) = v iatpol (X,t). (4) 

Eq. 2] is solved by a midpoint method which is straightforwardly integrated into the leapfrog scheme that advances 
the turbulent fields. Test calculations using Fourier interpolated 'exact' representations of turbulent velocity fields 
have shown that the chosen tricubic polynomial interpolation delivers sufficient precision with a mean relative error at 
512 3 resolution of ~ O(10 -3 ). In addition tricubic interpolation is numerically much cheaper than the nonlocal spline 
approach (cf. [2lj|). especially on computing architectures with distributed memory. The initial particle positions 
are forming tetrads that are spatially arranged to lie on a randomly deformed cubic super-grid with a maximum 
perturbation of 25% per super-grid cell. This configuration represents a compromise between statistical independence 
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TABLE I: Particle groups and respective initial pair separations Ao 



of particle dynamics and a space-filling particle distribution (cf. [ij, l22j]). In addition well-defined initial particle-pair 
separations A (cf. Table |J) are realized by the tetrad grouping. 

Distances are given in units of the Kolmogorov dissipation length 4 = (/^ 3 /e K ) 1 / 4 defined with the kinetic energy 
dissipation rate e K that is part of the total dissipation rate e = e K + e M = J v dV(fiu} 2 + ry(V x b) 2 ). In this work 
^max^d ~ 1-6. This fulfills the widely accepted resolution criterion introduced in and corresponds to a dissipative 
energy fall-off of about 4 decades in the dissipation range avoiding interpolation problems at grid-scales. Intervals of 
time are given in units of the large-eddy turnover time, X = ir/f^E) 1 / 2 , E = E K + E M = j v dV(v 2 + b 2 )/2 being 
the total energy, or in multiples of the Kolmogorov time, Td = (^/e 1 ^) 1 / 3 , as appropriate. 

The simulations are carried out using a resolution of 512 3 collocation points with aliasing errors being treated by 
spherical mode truncation [23| . Quasi-stationary turbulence is generated by a forcing which freezes all modes in the 
sphere k < kt = 2. The frozen modes which are taken from DNS decaying turbulence sustain the turbulence gently 
via nonlinear interactions. It has been checked that this way of driving does not introduce significant anisotropy by 
regarding direction-dependent Eulerian two-point statistics. 

Starting with a set of random fluctuations of v (and b in the MHD case) with zero mean the driven flows reach 
quasi-stationary states during which the total energy E shows fluctuations < 10% around unity and E M /E K as 2 
(MHD). In both simulations the total energy dissipation rate e is quasi-constant at about 0.24 with e M « 0.15 in the 
MHD case. The turbulent fields in the MHD system have negligible magnetic and cross helicity. The macroscopic 
Reynolds numbers are dimensionaily estimated using fj,, 77, e, E, and the kinetic energy E K as Re = (E K y/ 2 E 3 / 2 /(e//) 
(hydrodynamic) and Rm = RePr m (magnetic) with the magnetic Prandtl number Pr m = fi/rj set to unity. The 
respective numerical values of the parameters are /1 = 8 x 10~ 4 , Re w 5400 (NS) and /1 = r\ = 5 x 10~ 4 , Re « 5200 
(MHD). Cases with Pr m ^ 1 while interesting due to their importance in the context of turbulent dynamos (see e.g. 
24l 25j |) are beyond the scope of this paper and will be addressed in future work. 



After the runs have reached macroscopic quasi-equilibrium the trajectories of massless point particles marking the 
fluid are traced over 13.5X (NS) and 11T (MHD) corresponding to 500rd and 350rd, respectively. The initial particle 
positions are chosen in five groups of tetrads with different particle-pair separations (cf . Tabic fl} . 

In statistically isotropic turbulence single tracer particles are expected }26j to show a diffusive time dependence of 
the mean-square particle displacement {(X(t) — -X^o) 2 ) ~ t, Xq = -X"(0), for t 3> XL where XL is the autocorrelation 
time of the Lagrangian velocity V(t). Here, XL ~ I6rd (NS) and Xl ~ 15rd (MHD). If t <C XL ballistic scaling is 
predicted, ((X(t) - X ) 2 } - t 2 . 

In both simulations (cf. fig. [J) ballistic scaling can be identified up to about Xl. Diffusive behavior is observed 
for t > 50rd. At t > 70rd the particles have traveled about 5004, i.e. half the size of the simulation volume, and 
finite-size effects as well as the influence of the large-scale driving can be detected. The normalized turbulent diffusion 
coefficient Aurb(**) = / dr(V(t)V(t + T))/{V 2 (t)) with the avergage running over all trajectories shows for both 
systems in the interval < t* < 50rd a steep increase with a subsequent saturation at the asymptotic value Xl- 

It is found that with regard to turbulent single-particle diffusion the NS and the MHD system show no significant 
differences. The small offset of the MHD displacement curve compared to the NS simulation is explained by the lower 
level of kinetic energy in the MHD system which is not fully compensated by the applied normalization. An analytically 
predicted slowing down of diffusion (and dispersion) [2?| is not found here. The cited result is however based on the 
restricting assumption of a velocity field which is delta-correlated in time thereby neglecting the dynamically important 
adaptation of the velocity fluctuations to the magnetic field structure. 

The observed similarity of the curves in fig. [1] indicates that statistics of single-particle trajectories is not a 
proper instrument to study the structural differences in the velocity field of the NS and MHD systems caused by 
macroscopically isotropic magnetic field fluctuations (cf. 0, 0, [3(| for numerical simulations). 

In this respect a more instructive diagnostic is relative pair-separation or dispersion statistics where the separation, 
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FIG. 1: Evolution of normalized mean-square distance to initial position for turbulent single-particle diffusion in Navier-Stokes 
(black) and MHD turbulence. The dashed lines indicate ballistic scaling ~ t and diffusive behavior. 




FIG. 2: Evolution of normalized mean-square relative dispersion in Navier-Stokes (black) and MHD turbulence. The behavior 
for three different initial pair-separations (particle groups in Table[I} are shown, solid: group 1, dashed: group 3, dash-dotted: 
group 5. Straight lines denote ballistic ~ t 2 and diffusive scaling. 



A = X±(t) — X%(t), of two particles is considered 3l|. The pair-separation in the ballistic regime, t <C tj, where 
the particle velocities are finite and constant obeys the relation ((A(i) — A ) 2 ) ~ t 2 (see @ for experimental and 
[9l. fioj] for numerical results (all NS)) while for t 3> Tl diffusive scaling holds, ((A(t) — A ) 2 } ~ t since the dynamics 
of particles forming a pair are statistically independent in this case. 

The asymptotic limits can be identified in fig. [2] with differences between NS and MHD statistics showing up at 
intermediate times. For clarity, only data based on the particle groups 1, 3, and 5 (cf. Table HJ are shown, fn both 
systems the evolution of groups 2 and 4 is qualitatively similar to the behavior of group 3. As expected ballistic 
scaling ~ t 2 is visible in both simulations for short times up to about ra- Eventually, for t > 160rd an approach 
to the diffusive limit is seen rudimentally. The dispersion is accelerating during 10 < t < 100. However, neither 
the Batchelor law [Hj], ((A - A ) 2 ) = ll/3C 2 (eA ) 2/3 i 2 , (cf. (H) nor Richardson scaling H, (A 2 ) - et 3 , (cf. 
SEH) can be clearly identified since the simulations do not generate sufficiently large inertial ranges. In addition, 
the theoretical preconditions necessary for Batchelor and Richardson behavior are not satisfied in the simulations 
presented here. In particular, Batchelor dispersion requires Ao to lie in the inertial range 31( for recovering of the 
exact prefactor. Richardson scaling is also not expected since it would entail very large A in the inertial range [33| . 
In addition Richardson behaviour would imply an appoach of the pair-separation curves to one universal scaling law 
independent of Aq which is not observed here. 
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FIG. 3: Normalized and averaged pair-separation velocity, {Uu), in the direction of A using the same symbols as in fig. [5] 
Black lines: NS, grey lines: MHD 



While the evolution of NS and MHD dispersion is qualitatively similar the acceleration phase in the MHD system 
is significantly delayed compared to the NS case. The reason for this difference is found in the averaged separation 
velocity (U\\) in the direction of the separation vector A with U = A = C/yA/A + U± shown in fig. [3] in units 
of Vd = £d/ T d- The separation velocity for all particle groups except groups 5 (cf. Table Q} which have the largest 
initial separation, A , displays a continous increase before passing through a slowing-down phase for t > t<j. The 
beginning of the slowing-down marks the point of time at which the particles start to sense temporal fluctuations 
of the velocity field. The subsequent acceleration of dispersion can be understood by the increasing importance of 
sweeping by more coherent larger-scale eddies. The maximal separation velocity for all particle groups except groups 
5 is reached around 90rd (NS) and 150rd (MHD). There the mean pair separation is about half the extent of the 
periodic simulation volume and the particles start to approach each other again. The temporal shift of the MHD 
maxima compared to the NS curves as seen in fig. [3] is explained by the smaller kinetic energy of the MHD system 
and the observed stronger slowing-down of the average MHD pair-separation velocity. The separation velocity curves 
for the largest initial pair-separations (groups 5) do not display the slowing-down phase since these pairs probe only 
the largest spatial scales of the flow where the driving is governing turbulent dynamics, cf. 

The main difference between the NS and MHD cases, however, lies in the slowing-down phase which is much more 
pronounced in the MHD simulation. 

The reason is the well-known anisotropy of turbulent eddies with respect to the local magnetic field, see for example 
34 . 35 , d^, [I?], [H, . Since small-scale eddies are elongated in the local magnetic field direction MHD fluid elements 
are more likely to travel in similar directions oriented along the local magnetic field. The field-parallel velocity 
component is causing the effective particle pair-separation while the field-perpendicular components are associated 
with Alfvenic quasi-oscillations which do not contribute to the average particle displacement. Contrary to the ballistic 
regime with quasi-constant flow velocities for t <§; Td, the anisotropy of the fluctuating velocity field at later times has a 
constricting effect on turbulent dispersion compared to the NS case. Consequently, the fluid particles are preferentially 
traveling along the magnetic lines of force which significantly reduces dispersion in the field-perpendicular directions. 

This conjecture is supported by fig. [4] which shows probability distributions of the angle 7 = /(£/, B) for particle 
group I and points of time in the interval < t < 64rd introducing a rough proxy of the mean magnetic field at scale 
A, B = [b{X 1 ) + b{X 2 )}/2. 

For isotropic random velocities a sinusoidal distribution (thin solid line) would be expected. However, it is seen 
that even for comparably large times the distribution of the angle P{^) exhibits a clear deviation from this behavior 
favoring velocities aligned with the magnetic- field proxy B. The observed trend to sinusoidality with increasing time 
is due to (A) approaching the largest scales of the flow which leads to weakly correlated fluctuations in U and B. 
This trend is limited by the way of forcing chosen in this work. 

Apart from pair-dispersion, tracer dynamics display another interesting difference between Navier-Stokes and MHD 
turbulence. At smallest scales in the vicinity of the most singular dissipative structures the tracer trajectories differ 
significantly. In the NS simulation where the smallest flow structures are vortex filaments the fluid particles describe 
helical motions around the filaments (cf. also [l3|). In contrast, vorticity sheets typical for smallest-scales in MHD 
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FIG. 4: Probability distributions of the angle 7 between the relative velocity U = A of particles in group 1 (cf. Table [T]) and 
the proxy B = [b(Xi) + b(X2)]/2 at different times t (ra: long-dash, 4rd: dash-dot, 8rd: short-dash, 64Td: solid). The thin 
line denotes the distribution P(~y) = (7r/360) sin(7) indicating isotropic dispersion. 

turbulence lead to characteristic kinks in the tracer path. While the resulting trajectories are strongly different and 
characteristic for the respective turbulent system, these events occur too seldomly to have a measurable effect on the 
statistics of diffusion and dispersion regarded in this paper. 

In summary it was shown by comparison of direct numerical simulations of macroscopically isotropic Navier-Stokes 
(NS) and nonhelical magnetohydrodynamic (MHD) turbulence that the magnetic field in MHD turbulence slows down 
average pair-dispersion for intermediate times, < t < 10rd, compared to NS behavior. This effect is shown to be due 
to alignment of turbulent velocity and magnetic field fluctuations. Significant differences in turbulent single-particle 
diffusion could not be detected. Fluid particle trajectories in the vicinity of the strongly dissipative structures are 
characteristically different although these events have a negligible influence on the statistics investigated in this work. 

The authors would like to thank Holger Homann and Rainer Grauer (Ruhr-Universitat Bochum) for stimulating 
discussions and gratefully acknowledge support by Lorenz Kramer and Walter Zimmermann (Universitat Bayreuth). 
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